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ABSTRACT 

Recent tentative findings of non-Gaussian structure in the COBE-DMR dataset 
have triggered renewed attention to candidate models from which such intrin- 
sic signature could arise. In the framework of slow roll inflation with built-in 
non linearities in the inflaton field evolution we present expressions for both 
the cosmic microwave background (CMB) skewness and the full angular bis- 
pectrum Ci-^i^is in terms of the slow roll parameters. We use an estimator for 
the angular bispectrum recently proposed in the literature and calculate its 
variance for an arbitrary ii multipole combination. We stress that a real detec- 
tion of non-Gaussianity in the CMB would imply that an important compo- 
nent of the anisotropics arises from processes other than primordial quantum 
fluctuations. We further investigate the behavior of the signal-to- (theoretical) 
noise ratio and demonstrate for generic inflationary models that it decreases 
in the limited range of small-£'s considered for increasing multipole i while 
the opposite applies for the standard 's. 

Key words: cosmic microwave background - methods: analytical - cosmol- 
ogy: theory - large scale structures of Universe - early Universe. 



1 INTRODUCTION 

The theory of inflation provides an elegant means to solve the usual problems (horizon 
and flatness problems) of the standard model of Cosmology (Guth 1981). It consists in 
assuming that a phase of accelerated expansion took place in the very early Universe at 
the GUT energy scale. In the most simple models of inflation, this phase of accelerated 
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expansion is driven by a scalar field. The physical origin of this field is still an open problem 
and Physics at GUT scale could be much more complicated than assumed in these simple 
models. However, it should be stressed that the very concept of inflation lies in the fact that 
the second cosmic time derivative of the scale factor was positive in the early Universe and is, 
in this sense, independent of any model-building provided that the effective inflaton potential 
is flat enough to allow an inflationary expansion of at least ~ 70e-folds. Therefore slow roll 
inflation can be viewed as a generic framework which permits to implement concretely the 
concept of inflation and allows to perform simple analytical calculations. 

The beauty of the inflationary scenario is that, combined with Quantum Mechanics, 
it also provides a natural explanation of the origin of the large scale structures and of the 
cosmic microwave background (CMB) radiation anisotropics observed in our Universe (Guth 
&Pi 1982; Starobinsky 1982; Hawking 1982; Bardeen, Steinhardt & Turner 1983). In this 
explanation, the quantum character of the inflaton field plays a crucial role since the seeds 
of these perturbations are the unavoidable quantum fluctuations present at the beginning 
of the inflationary epoch. Then these fluctuations are parametrically amplifled during the 
accelerated phase of expansion (Grishchuk 1974). Therefore the properties of the initial 
spectrum of perturbations depend on the initial quantum state in which the fluctuations 
were placed and on the behaviour of the scale factor during inflation. Observationally, we 
have access to the initial spectrum when one looks at the CMB anisotropy multipolc moments 
corresponding to the largest angular scales on the celestial sphere. Indeed these multipoles 
are dominated by modes whose wavelengths are comparable to the size of the horizon today. 
This means that after their creation these modes spent most of their time outside the Hubble 
radius and as a consequence were not contaminated by astrophysical processes: in a certain 
sense, they can be viewed as a pure relic of the very early Universe. 

Among the many features of the perturbations, the statistical properties are certainly of 
a big importance. In the theory of cosmological perturbations of quantum mechanical origin, 
it is assumed that the initial state is the vacuum. This seems to be the most natural choice 
although it was already noticed that it could be difficult to understand why the fluctuations 
of a fleld which is initially out of equilibrium would be placed in this state (Unruh 1998). 
It was recently shown (Martin, Riazuelo & Sakellariadou 1999) that if one tries to start the 
evolution from a non vacuum initial state, then observations require that this state be close 
to the vacuum. This seems to indicate that the vacuum is indeed a reasonable choice. Since 
each mode of the perturbations can be viewed as an oscillator, one immediately reaches the 
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Cosmic microwave background bispectrum and slow roll inflation 3 

conclusion that the corresponding statistical properties must be Gaussian (recall that the 
ground-state wavefunction of an harmonic oscillator is a Gaussian function). This constitutes 
an important and generic prediction of inflation. 

Another source of cosmological relevant density inhomogeneities arises in models with 
topological defects of the vacuum, like cosmic strings and textures. These would leave dif- 
ferent imprints on the CMB both at recombination and later, during the photon travel from 
the surface of last scattering to the present, on various angular scales (Allen et al. 1997; 
Magueijo et al. 1996; Battye et al. 1998; Contaldi et al. 1999; Pogosian & Vachaspati 
1999; Durrer et al. 1996; Pen et al. 1997). There is the hope that future balloon-borne 
(e.g., MAXIMA, Lee et al. 1999; BOOMERanG, Lange et al. 1999) and satellite (MAP 
and Planck surveyor, Bersanelli et al. 1996) missions will allow a clean distinction among 
these different classes of models by mapping the CMB with unprecedent precision. 

A number of authors (e.g., Smoot et al. 1994; Torres et al. 1995; Hinshaw et al. 1994; 
Hinshaw et al. 1995; Kogut et al. 1996) have analysed the COBE-DMR sky maps with a 
variety of test (like three-point statistics, genus, and extrema correlation function) and found 
perfect agreement with a Gaussian distribution. Recently however, three groups (Ferreira 
et al. 1998; Pando et al. 1998; Bromley & Tegmark 1999) have analysed the COBE- 
DMR four-year dataset and reported detections of non-Gaussianity casting doubts on these 
early findings. Banday et al. 1999 have further analyzed the same data, finding that the 
non-Gaussian signal is driven by the 53 GHz sky maps. They concluded that this frequency 
dependence strongly indicates that the signal is not of primordial origin. Despite all this, it 
seems now that due to limited signal-to-noise, sky coverage, and uncertainty in foreground 
substraction, present day experiments cannot conclusively exclude non-Gaussianity to a 
satisfactory confidence level. 

The above remarks regarding the Gaussian character of the primordial perturbations have 
been established within the framework of the linear theory of cosmological pertubations. It is 
clear that generic higher order, in particular quadratic, terms are present and will produce a 
non vanishing signal even for inflationary models (Linde & Mukhanov 1997). The predictions 
for the three-point correlation function on large angular scales due to nonlinearities in the 
inflaton evolution were considered in the past (Falk et al. 1993; Gangui et al. 1994; Gangui 
1994). 

The post recombination Rees-Sciama effect, due to the mildly non-linear evolution of the 
perturbations also contributes to the signal (Luo & Schramm 1993; MoUerach et al. 1995; 
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Munshi et al. 1995). Evolving networks of topological defects continuously seed perturba- 
tions on the CMB that, by the very nature of the sources, are predicted to be highly non- 
Gaussian (Bouchet et al. 1988; AveUno et al. 1999; Gangui & Perivolaropoulos 1995; Gan- 
gui & MoUerach 1996). Further secondary effects contribute to produce non-Gaussianities 
at smaller scales (Aghanim & Forni 1999) and would be characterized by detectable cor- 
relations between gravitational lensing and Sunyaev-Zel'dovich maps (Spergel & Goldberg 
1999; Goldberg & Spergel 1999). 

We here present a general discussion of non-Gaussian features arising in the framework of 
slow roll inflation. Our article borrows some definitions and formulas for the CMB three-point 
correlation function (specially in section 2) from (Gangui et al 1994). Our main aim here 
is to present explicitly the derivations of the non-Gaussian estimators as a function of the 
inflationary slow roll parameters in the Legendre space, namely the full angular bispectrum 
Cf^£2^3. Then, in the third section, we present for the first time the analytical expression for 
the variance of an estimator recently proposed for the bispectrum (Ferreira et al. 1998) in 
the mildly non-Gaussian approximation. This allows us, in the last section, to compare the 
behaviour of both quantities for various multipoles and conclude with the by now established 
result supporting the view that the recently observed non-Gaussianity cannot be explained 
in the framework of slow roll inflation. 



2 NON-GAUSSIAN SIGNAL IN REAL AND LEGENDRE SPACES 

In the framework of the theory of cosmological perturbations of quantum mechanical origin, 
AT/T{x,^) is an operator. The corresponding statistical properties arc then calculated 
by "sandwiching" this operator (or a combination of these operators) with the quantum 
state 1^) in which the quantum perturbations are placed. However, it has been shown in 
(Grishchuk & Martin 1997) that one can think to AT'/T(x, 7) as a classical stochastic 
process. This stochastic process can be expanded in spherical harmonics 
AT 

— (f,7) = E«r(^)>v,yr(7), (1) 

where Wi represents the window function of the particular experiment. The coefficients 
a™ (5) are random variables and are in principle different for different observers at positions 
X. The statistical properties of AT/T{x, 7) are completely specified if the probability density 
function (pdf) of the a™(x)'s is known. Choosing the initial state to be the vacuum, i.e. 
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|\E') = |0), and considering only linear terms is equivalent to saying that the pdf of the 
a^(x)'s is a Gaussian distribution (Grishchuk & Martin 1997). This means that 



where brackets (■) stands for an average over the ensemble of possible universes in the 
sense explained above. The variance is rotationally invariant, i.e. depends only on i, signal- 
ing statistical isotropy. For Gaussian fields the previous equations are sufficient since this 
kind of fields are completely characterized by giving their two-point correlation function or, 
equivalently, their (angular) power spectrum. 

We now take into account the non linearities. This means that the pdf of the a™(x)'s 
is no longer a Gaussian distribution. The two first moments will still be given by Eq. 
and the difference will show up at the level of the third order moment. Predictions from 
different models usually come as expressions for the ensemble average which 
can be written in full analogy with Eq. in terms of the angular bispectrum Q^^j^g as 
follows: 



where now the proportionality factor is a Wigner 3-j symbol. This is non-zero only if the 
indices (i = 1,2,3) fulfill the relations: \ij — ik\ < < + and mi -|- m2 + 

m3 = 0. There is an additional "selection rule" in this equation that arises from demand- 
ing that (AT(7i)AT(72)AT(73)) be invariant under spatial inversions. One then obtains 



For a Gaussian model we clearly have Q^^a^s — whereas the non linear evolution of the 
perturbations will induce a Qi^a^s 7^ 0- Therefore in order to probe the Gaussian character 
of the stochastic process AT/T(x, 7), it is certainly convenient to study quantities related 
to the third moment. A priori a large choice is allowed. Here below we will concentrate on 
the CMB collapsed three-point correlation function, defined in real space. The skewness will 
just be the particular case of the collapsed function at zero lag. 

The three-point correlation function for points at three arbitrary angular separations a, 
(3 and 7 is given by the average product of temperature fiuctuations in all possible three 
directions with those angular separations among them (Gangui et al. 1994). The collapsed 
case corresponds to the choice a = j3 and 7 = and reads 



(a-(f)) = 0, (a™^(f)a,7*(x)) = C,,5,,,,5, 



(2) 




(3) 



(a^^a^^^a?^) = for + £2 + 4 = odd (Luo 1994). 




(4) 
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For a = 0, we recover the well-known expression for the skewness, C^lO). Using the spherical 
harmonics expansion (|l]) the last equation can be rewritten as: 

£1/2 /a, 1^1, 1^2, m3 

where we have defined the coefficients Ti-^l^l^^ 

HT,}7r = / d^^K' (7)>^r (7) , (6) 

which has a simple expression in terms of Wigner 3-j symbols (Messiah 1976): 



^mimama _ , / (2^1 + l)(2f2 + 1)(2£3 + 1) fe, es\ /£i h h \ (7\ 
'^hl2t3 ~ \l A„ lo Mmi ma ma i • \' ) 



An 

We see here that the condition ii + £2 + (^3 = even is enforced by the presence of the first 
Wigner 3-j symbol. Substitution of Eq. (|^) together with Eq. (|^) into Eq. where we have 
taken the ensemble average, yields for the mean collapsed three-point correlation function 

We see from this that all terms in the sum satisfying ii + £2 + £3 = even and the triangle 
inequalities but otherwise arbitrary will contribute to the value of the collapsed three-point 
function and hence also to the skewness. In general a complete probe of the three-point 
function will require the knowledge of all the coefficients ^^£2^3 cind not just the "diagonal" 
ones, Cm. 

To estimate the amplitude of the non-Gaussian character of the fiuctuations one usually 
considers the "dimensionless" skewness Si = {Cs{0)) / {C2{0))^^'^ ■ Alternatively, if we want 
our results to be independent of the normalisation, we may also define the ratio ^2 = 
(C3(0))/(C2(0)) . Both quantities will be known once the bispectrum has been calculated. 
Therefore our aim is now to compute ^^£2^3 in the context of slow roll infiation, i.e. in terms 
of the infiationary slow roll parameters. 

In the framework of the stochastic approach to infiation (Starobinski 1986; Goncharov, 
Linde & Mukhanov 1987), the calculations reported in (Gangui et al 1994) valid for models 
satisfying slow roll dynamics yield 

where in general one requires that the infiaton potential V{(f)) be a smooth function of its 
argument, which translates into requiring well defined values for the steepness of the potential 
X = nipiV'/V (here ' = d/d0 and mpi is Planck mass) and its derivatives throughout the 
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range of relevant scales (Turner 1993). Another way of expressing this result is in terms of 
the standard slow roll parameters (Dodelson et al. 1997; see also Stewart & Lyth 1993; 
Liddle et al. 1994) 



^2 /v'\^ 





'V" 




?)1 
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-a 





(10) 



For the very large scales we focus on here and for standard chaotic initial conditions in 
the inflaton field, both parameters satisfy e, <^ 1. In terms of these, a comparison of Eqs. 
(I) and (H) leads to: 
5 



Ce,e,e, = ^(3e - 2r])^{2i, + + 1)(2£3 + ^tA. + C^U, + C.-M ('o^ %) ■ (11) 

This equation shows explicitly the full angular bispectrum arising from generic inflationary 
models in terms of the slow roll parameters. To be specific in the sequel we will take a 
potential V^(0) oc 0^, with p > 1. Then, after evaluation of the slow roll parameters at Hubble 
radius crossing for the relevant very large scales one has (e, r/) = {pj (p+200), (p— 2) / (p+200)) 
corresponding to a scalar spectral index n = 1 — (2p + 4) / (p+200). The calculations reported 
in the following sections will be performed for a quadratic potential, namely, p = 2. 

With this result in mind, we can now re-express the three-point correlation function 
and skewness. We first need to evaluate the multipole moments Ci. For large scales we have 
Pq>{k) oc k^~^ with n corresponding to the primordial index of density fluctuations (e.g., 
n = 1 is the Harrison-Zel'dovich, scale invariant case) in which case (Bond & Efstathiou 

1987; Fabbri, Lucchin & Matarrese 1987) 
^ ^ T{^ + n/2 - l/2)r(9/2 - n/2) _^ ~ 

- T{1 + 5/2 - n/2)r(3/2 + n/2) = ^^^^ 

with C2 related to the quadrupole power spectrum normalization Qrms-ps = ^o(5C2/47r)^/^. 
Sometimes it turns out to be convenient to factorise the quadrupole amphtude out by us- 
ing Cg. In general, the quadrupole depends on the spectral index n. Thus, normalization 
analyses of datasets yield the estimate of the pair {n, Qrms-ps)- For example, the maximum 
likelihood analysis of the COBE-DMR dataset performed in (Bunn & White 1997) yields 
(n, Qrms-ps) = (1-2, 16.2^K) while their best fit scale invariant normalization that we will 
use for the numerics in the next section is Qrms-ps = ^8.7fiK (same as in Gorski et al. 
1996). 

It is now easy to obtain Si and ^2 in terms of the slow roll parameters. They are given 

by 
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J r(3-n)r(3/2 + n/2) 
= mi,i — ^^' - i p(2 - n/2)r(9/2 - n/2) ) ^^^^ 

and 

52 = 15(3e-2r7)J2(n), (14) 

where the normahzation dependence (in Si) is made exphcit by defining the spectral index- 
dependent geometrical factor 



Tain) 



|E^,/,.^3(2^i + l)(2^2+l)(2£3+l)[C£.C£.+C£.Q3+^4^^i]>V^i>V^2>V^3-^4^2^3 



(15) 

E^(2/ + i)Qwi; 

The exponent q in the denominator takes values 3/2 and 2 for Si and respectively. 
The coefficients Ti^^^i.^ = (47r)~^/ rffi^ / (ifiyP^^(7 ■ 7')P£2(7 ■ 7')P^3(7 ■ 7') may be suitably 
expressed in terms of products of factorials of £1, ^2 and £3, using standard relations for 
Wigner 3-j symbols: in fact we have J^Mm =(00™)- COBE-DMR window function, 

the numerical factors Tq{n) in Eq. ( plSf ) are of order one for all interesting values of the 
primordial scalar spectral index. Eqs. (|l^) to (|15D were already presented in a different 
form in (Gangui et al 1994). Particular cases of these equations have also been displayed in 
(Kamionkowski & Kosowsky 1999; Verde et al. 1999). 



3 BISPECTRUM ESTIMATOR AND ITS VARIANCE 

When one particular mechanism for the generation of CMB non-Gaussian features is speci- 
fied, it is a direct procedure to compute the analytical angular bispectrum. One such example 
was shown in the previous section in the case of slow roll infiation. However, when dealing 
with just one realization of a stochastic process, as is the case for the CMB, all computed 
quantities come with theoretical error bars (Scaramella & Vittorio 1991; Srednicki 1993). 
Even though we can analytically compute mean values, when an actual observation is made 
there is a non-vanishing probability that it will fall within a value ±0" apart from the mean. 
This problem has been dubbed "cosmic variance". To deal with it, one has to introduce an 
estimator E of the quantity e we seek, i.e. a random variable such that {E) = e. In this case 
the estimator is said to be unbiased. Then one should compute the variance of the estimator, 
a^, and try to minimize it. If it turns out that = then we can find e with the help of one 
realization only (in fact because each realization gives e). In general, we have 7^ 0, and 
one can show that this is linked to the fact that a stochastic process cannot be ergodic on a 
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(celestial) sphere (Grishchuk & Martin 1997). In that case will express the unavoidable 
error made when one estimates the mean of a stochastic process from one realization. 

For the standard angular spectrum Ce the best unbiased estimator is (Grishchuk & Martin 
1997; Tegmark 1997) 

h = ^ E (16) 

and it is easy to check that it is unbiased, namely (fg) = Ce. Its variance, the smallest one 
amongst all possible estimators variances, is given by 



(^f = \Mr—Ci- (17) 

It is clear that such an optimal strategy should be followed for the bispectrum as well (and 
for the higher order moments). The following expression 

f — \ " h t-z \ „rax mi^mz /-lo^ 

Je.xe.2h - 1^ \ra^ ma msj '^£1 ^l2 ^Iz \^^) 

is an unbiased estimator of the bispectrum C^^^j^g since we easily check that (/^i^a^s) = C^i^a^a- 
Its variance (squared) 



a 



2 _ /^2 \ /} \2 



(/<?1^2^3) " (/^ife^s) (19) 



will give us a first indication of the theoretical uncertainties we have to deal with. However, 
it should be clear that at this level nothing tells us that the one given in Eq. ([18|) is the best 
unbiased estimator for the bispectrum. As a consequence, working with it might be not the 
best choice and its variance might well be not the smallest one. Finding the best estimator 
for the bispectrum is not a trivial task and is presently under investigation (Martin & Gangui 
1999). Foregrounds, detector noise, sample variance in the cut sky are among the additional 
issues that need be mastered before claiming a real non-Gaussian detection. 

Recently, similar analyses for the computation of the variance for the estimator of 
(a^^a^^a^^) were presented (Luo 1994; Heavens 1998). If compared with the analysis 
of (Luo 1994) note that we are not estimating a^^^a^^^a^J^ but the statistically isotropic 
combination of Eq. (^) and there the presence of the 3-j symbol makes the whole differ- 
ence. Positive detection of intrinsic non-Gaussianity in the COBE-DMR four-year dataset 
was recently suggested (Ferreira et al. 1998; Pando et al. 1998; Bromley & Tegmark 1999). 
In particular, Ferreira and collaborators, in the attempt to unveil an eventually obscured 
non-Gaussian signal in real space, worked in the Legendre space and made use of an esti- 
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mator in the lines of Eq. ( jlSD above but with £i = £2 = ^3 = ^- They also normalized it by 
dividing fm by the estimator of [/^ in the notation of Eq. ([T6|)] to the power 3/2. 

We expect departures from Gaussianity to be weak and hence neglect the contribution 
of {feii'zk)"^ Eq. (|19|). In this mildly non-Gaussian approximation and after some straight 
algebra we obtain 

■'«1^2^3 

where we demanded ^1 + ^2+^3 = even (otherwise Q^^j^g = 0) and 7^ 0, what considerably 
simplified the resulting expression. We note in passing that for the above computation it is 
useful to recall the identity (Mollerach et al 1995) 

E (-1)™ U L ''0) = i-iYV2lTi 4,0. (21) 

m=-e 

We are now in a position to compare the signal, i.e. the bispectrum given in Eq. ([Tip, with 
the theoretical noise characterized by (|20|). 



4 DISCUSSION AND CONCLUSIONS 

In the previous sections we have computed both the expression for the angular bispectrum, 
as obtained generically in the framework of slow roll inflation whenever one goes beyond 
the linear order, and the variance associated with an unbiased estimator, assuming a mildly 
non-Gaussian process. We can now compare these results for an arbitrary configuration of 
ii multipoles. As a representative example, and given the fact that this was actually the 
case considered in the literature, we consider ii = £2 = £3 = i = even. We then have the 
bispectrum 

Cm = ^(3e - 2v){2£ + 1)^/^ C] ^ , (22) 
while the variance is now given by 

af^^^ = V6Cf. (23) 

We show the relative amplitudes in Fig.l. The plot allows us to judge how plausible it is 
for generic one-field inflationary models to reproduce any possible non-Gaussian structure 
found on large angular scales, in particular on the COBE-DMR dataset. Single different 
values (correlations) for the indices £±,£2, £3 (satisfying £1 + £2 + £3 = even and the triangle 
inequalities) can be tried with similar result. Leaving aside for the time being the possibility 
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Figure 1. Normalised angular bispectrum Cm as predicted by a generic slow roll inflation model, in units of C|, with C2 = 
1.18 X 10"^" related to the quadrupole power spectrum normalization Qrms— PS = l8.7fiK, as a function of the multipole index 
£ for all even values up to £ = 20 (white dots in main plot). Grey band corresponds to the normalized variance f^f^^^ (also in 
units of C|) associated to the estimator of Eq. (|l^). In the inset we zoom up l'^{£ + l)-^Cut in the same units, which permits 
to see the alternating sign of the normalised bispectrum and its actual smooth increase in amplitude with increasing i. 



of foreground contamination and assuming any non-Gaussian signal is intrinsic to the CMB, 
we should conclude that the presently considered class of models cannot explain it. 

In order to be more specific, we now turn to the study of the "signal-to-noise" ratio 
defined by the following expressions 

\ Cf ^ S \ Ce 



(24) 

fi \^ » / 3,^ f^^^ 

for the angular spectrum and bispectrum respectively. For the first one [neglecting speci- 
ficities of the particular experiment, detector sensitivity, pixelization, etc. (Knox 1995)] the 
following well-known behaviour is found 

S_\ _ / 2f + 1 
N)21~ V 2 



(25) 



This means that the signal emerges from the noise while going towards big values of £, 
accounting for the fact that the cosmic variance is important at large scales only. For the 
bispectrum, one has 
S_\ _ 15 



(3e-2r7)(2£+ l)2/2(- 



1/2 (t t t 

oy ; 



(26) 



2V67r 

and we see that, contrary to the previous case, the signal-to-noise ratio depends on Ci and 
on the slow roll parameters. The behaviour of |(S'/A^)3^£| is displayed in Fig. ^ where we see 
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Figure 2. Absolute value of the signal-to-noise ratio for the bispectrum {S/N)^ i versus the multipole index £. 

that it diminishes in absolute value with increasing i. This behaviour is easily understood 
once we look at the hierarchical form of ^^^3^3 and compare with the dependence of the 
variance. Again, even restricting ourselves to small £'s, we see that the presence of the Wigner 
3-j symbol is the responsible for this particular behaviour. Hence, approaching the largest 
multipoles in the COBE-DMR data set (in particular i = 16) the situation gets worse. In 
fact, from Eq. ( p6|) we roughly have Ci oc + 1)~^ and i q q) almost constant with 
i (see Fig. |]), and then \{S/N)s^i\ oc while = |(5'/A^)2/| oc in the same range of 

validity. Note that the above analysis should be supplemented by a similar one wherein the 
behaviour of the bispectrum coming from the post recombination Rees-Sciama effect is also 
considered. However, given the smallness of the non-Gaussian signal this should not modify 
our conclusions very much. 

To conclude, let us emphasize again the main results obtained in this article. We have 
stressed that a real detection of non-Gaussianity in the CMB would imply that an important 
component of the anisotropics arises from processes other than quantum fluctuations during 
an early inflationary epoch. This notwithstanding, inflation predicts an actual generic form 
for the bispectrum, Eq.(|ll|), and we here showed it explicitly in terms of the slow roll 
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Figure 3. Absolute value of the product £ q for different £'s. Note tliat the product vanish for all £ = odd. 

parameters. We also computed for the first time the variance of one candidate estimator for 
the bispectrum often employed in the literature and showed that the signal is drowned in 
it. Contrary to the standard spectrum case, one cannot hope to palliate at least somewhat 
this problem by going to higher values of i (always within the small-£ region) since the 
signal-to-noise ratio decreases with i like However, it should also be stressed that this 

conclusion might well be weaken by the finding of the actual best unbiased estimator. We 
hope to address this question elsewhere. 

Note added: After the submission of this paper, a preprint by Wang and Kamionkowski 
( |astro-ph/9907431j ) appeared in which similar conclusions were reached. 
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